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We present finite temperature (T) extension of the QED3 theory of underdoped cuprates. The 
theory describes nodal quasiparticles whose interactions with quantum proliferated hc/2e vortex- 
antivortex pairs are represented by an emergent U(l) gauge field. Finite T introduces a scale beyond 
which the spatial fluctuations of vorticity are suppressed. As a result, the spin susceptibility of the 
pseudogap state is bounded by T 2 at low T and crosses over to ~ T at higher T, while the low-T 
specific heat scales as T 2 , reflecting the thermodynamics of QED3. The Wilson ratio vanishes as 
T — » 0; the pseudogap state is a "thermal (semi-)metal" but a "spin/charge dielectric." This non- 
Fermi liquid behavior originates from two general principles: spin correlations induced by "gauge" 
interactions of quasiparticles and fluctuating vortices and the "relativistic" scaling of the T — 
fixed point. 



Recent experiments 0- IH Q provide support for the 
view that the pseudogap state of copper oxides repre- 
sents a phase-disordered superconductor This view 
is central to the QED3 theory of underdoped cuprates 
whose degrees of freedom, Bogoliubov-deGennes (BdG) 
quasiparticles and fluctuating hc/2e vortex- antivortex 
pairs, and their mutual interactions are argued to capture 
the effective low energy physics of a d-wave superconduc- 
tor in a doped Mott insulator. Within the pseudogap 
phase of high-T c cuprates this QED3 theory assumes the 
role played by the Fermi liquid theory in conventional 
metals and superconductors. The theory possesses three 
major dynamical symmetries: relativistic and gauge in- 
variance and chiral symmetry, all three emergent in na- 
ture 0|. Irrespective of whether the symmetric phase of 
QED3 is the true T — ground state of the system, these 
symmetries conspire to make its peculiar brand of quan- 
tum criticality a strongly attractive basin of influence on 
physical properties of the pseudogap state. 

In this Letter we focus on the subset of such properties 
which arc intrinsically important features of a condensed 
matter system: thermodynamics and spin response of 
the pseudogap state. Our primary aim is to stimulate 
experimental activity by providing explicit and testable 
theoretical predictions stemming from the theory of Rcf. 

The summary of our results is as follows: first, in 
order to derive thermodynamics and spin susceptibility 
from QED3 theory we generalize its form to finite T. 
This is a matter of some subtlety since, the moment 
T ^ 0, the theory loses its fictitious "relativistic invari- 
ance". Second, we show that the theory predicts a fi- 
nite T scaling form for thermodynamic quantities in the 
pseudogap state and propose that this form be tested 
experimentally. Third, we determine the leading T 
scaling and demonstrate that the deviations from "rel- 
ativistic invariance" are actually irrelevant for T much 
less than the pseudogap temperature T*, in the sense 
that the leading order T ^ scaling of thermodynamic 



functions remains that of the finite-T symmetric QED3 
7J. These deviations from "relativistic invariance" do, 
however, affect higher order terms. Finally, we evaluate 
the uniform magnetic spin susceptibility x an d show that 
it is bounded by T 2 at low T but crosses over to ~ T at 
higher T, closer to T*. Consequently, the Wilson ratio 
R = \T/c v vanishes as T — > 0: the QED3 theory im- 
plies the non-Fermi liquid nature of the pseudogap state 
in cuprates. Such a state is a thermal (semi-)metal but 
a spin and charge dielectric. @ Our results are thus sug- 
gestive of the breakdown of Wiedemann- Franz law in the 
pseudogap state. 

The spin susceptibility of a d-wave superconductor 
vanishes linearly with temperature. A way to understand 
this result is to notice that the spin part of the ground 
state wavefunction, being a spin singlet, remains unper- 
turbed by the application of a weak uniform magnetic 
field. However, the excited quasiparticle states are not 
in general spin singlets and therefore contribute to the 
finite temperature susceptibility. Because their density 
of states is linear at low energies, at finite temperature 
the number of quasiparticles that are excited is ~ ksT, 
each contributing a constant to the Pauli-like uniform 
spin susceptibility x- Thus x ~ T. 

When the superconducting phase order is destroyed 
by proliferation of unbound quantum vortex-antivortcx 
pairs, the low-energy quasiparticles are strongly inter- 
acting. The interaction originates from the fact that it 
is the spin singlet pairs that acquire one unit of angular 
momentum in their center of the mass coordinate, car- 
ried by an hc/2e vortex. This translates into topological 
frustration in the propagation of BdG "spinon" excita- 
tions. As a result, non-trivial spin correlations persist in 
the excited states of the phase-disordered ci-wave super- 
conductor. At low temperature these correlations can be 
described by an emergent U(I) gauge field 0] and lead to 
suppression of x relative to its value for non-interacting 
BdG quasiparticles. We shall argue below that x T 2 . 
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Similarly, in a ci-wave superconductor, linear density 
of the quasiparticle states translates into a T 2 depen- 
dence of the low-T specific heat. When the interactions 
between quasiparticlcs are included the spectral weight 
is transferred to multi-particle states. Within QED 3 
theory, however, the strongly interacting 1R (infra-red) 
fixed point possesses emergent "relativistic" invariancc 
and the dynamical critical exponent z = 1. Further- 
more, the effective quantum action for vortices, deep 
in the phase-disordered pseudogap state, introduces an 
additional lengthscale, the superconducting correlation 
length £ T ± (labels r and _L stand for time- and space- 
like, respectively). At T = this scale serves as a short 
distance cutoff of the theory and is generically doping (x) 
dependent. We then argue that under rather general cir- 
cumstances the low-T (T « r) electronic specific heat 
scales as T 2 while the free energy goes as T 3 . We now 
proceed to substantiate the above claims. 

Vorticity fluctuations: 
Deep in the phase disordered pseudogap state the fluc- 
tuations in the vorticity 3- vector = e^ v \d v a\ are de- 
scribed by the following "bare" Lagrangian: 
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KT, K T T I b 2 , (1) 



where cu n is Matsubara frequency, K,K T are related to 
the finite superconducting correlation length of the pseu- 
dogap state as K T cx £ 2 C /£ T , K oc £ r |5j and we have 
also set vp — va = 1 for simplicity. f Tj x are gen- 
eral scaling functions describing how ^[a^] is modified 
from its "rclativistically invariant" T — form as the 
temperature is turned on. They involve only the ther- 
mal length ~ vp/T and K,K T and satisfy the condition 
/ r (0, 0,0,0) = /x(0, 0,0,0) = 1. Physically, the modifi- 
cations embodied in f T ± are due to changes in the pat- 
tern of vortex-antivortex fluctuations induced by finite 
T. The explicit expressions for f Ti ± depend on the de- 
tails of a particular model for phase disorder within the 
pseudogap state - as we emphasize below, however, our 
results are either completely insensitive to such details or 
reflect only the most general features of / r ,j_. 

To handle the intrinsic space-time anisotropy, it is con- 
venient to introduce two tensors 



x A^WrA k 2 ( x u n k v 



£> M „ — [ 5ij ^ 2 ) > 



and rewrite the gauge field action as 



(2) 



In the above equations fc M = {cj n , k) i.e. fc 2 = w, 2 + k 2 . 
It is straightforward to show that 

U° A = K T f T (~, £-,KT, K t t\ (k 2 + ljI), 



K° B = K T f T (^, —, KT, K T T ) ujI 
+ Kf/^,—,KT,K T T)k 2 

\ k Wn 



(4) 



The gauge field a M couples minimally to TV Dirac 
spinors representing nodal BCS quasiparticles (N = 2 
for a single CuC>2 layer). Consequently, the resulting La- 
grangian reads 



C = tp {il^d^ + 7 M a M ) -0 + £ [a M ] 



(5) 



where the summation over iV fermion flavors is under- 
stood. The integration over Berry gauge field repro- 
duces the interaction among quasiparticles arising from 
the topological frustration referred to earlier. 

Specific heat and QED3 scaling of thermodynamics: 
The only lengthscalcs that appear in the thermodynamics 
are the thermal length ~ Vp/T, and the superconducting 
correlation lengths K, K T . At T = 0, the two correlation 
lengths K, K T enter only as effective short distance cut- 
offs of the theory since the electronic action is controlled 
by the IR fixed point of QED 3 . These observations al- 
low us to write down the general scaling form for the 
doping-dependent free energy of the pseudogap state 



T(T;x) = -^ N (K T (x)T,K(x)T) 

Vp 



(6) 



where ^N(x,y) is the thermodynamic scaling function 
and £o(x) is the reference ground state energy. Note 
that we still maintain vp — va for simplicity and will 
discuss Dirac cone anisotropy shortly. In the above scal- 
ing form K T is directly related to the T — ► finite su- 
perconducting correlation length of the pseudogap state, 
£sc(x) 0- The ratio K T /K describes the anisotropy be- 
tween time-like and space-like vortex fluctuations and is 
also a function of doping x. The scaling expressions for 
other thermodynamic functions can be derived from (JJjJ 
by taking appropriate temperature derivatives. 

We are interested in the low temperature regime (T <C 
T*) in which the thermal length vp/T is much longer 
than K T and K (vf/T 3> £ sc (x)) or, equivalently, in 
the limit $at(x — > 0,y — * 0). This is just the limit in 
which the free energy T © approaches the free energy 
of the finite temperature QED3 and the precise form of 
l[TJl becomes unimportant as long as f T> ± remain finite as 
T — » 0. Consequently, $>n(%, y) should be regular at x = 
y = and takes on the value Cqed which is a universal 
pos itive numerical constant within QED3 HE3 (seeRef. 
[ill IT2I for discussion of related issues). Therefore, we 
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finally find that, at T < T*, T = -C QE uT 3 /v 2 F , S = 
3CqedT 2 /vp and c v = 6CqedT 2 jv\. 

The physical origin of the above form for the low tem- 
perature scaling can be illustrated by the expression for 
entropy consisting of the free fermion contribution and 
the part generated solely throughgauge field (vortex fluc- 
tuation) mediated interactions |l0|: 

S = NS° f + S mt = NC f T 2 - -^AF, (7) 
where, for N — 2, Cf — 9£(3) /8ttvfva and 



AF— 



X=A,B 



d 2 k 

(2n) 2 



tan 



' du w . 

^ COth( 2T )X 



Zmll£ re * 



TZen^ ret 



Note that we have traded the sum over Matsubara fre- 
quencies {tu n } for an integral over real frequency o;.|l3| 

The retarded fermion polarization function H^ ret 
can be written as NTVx ( [k| /T, uo /T) while the 
bare gauge field polarization n^ e * has the form 
KT 2 V\ (\k\/T,uj/T, KT, K T T) . By rescaling momenta 
and frequencies appearing in Eq. JSJ, we see that, 
in the low temperature limit KT <C 1 (T « T* or 
vf/T 3> £sc(a:)), the bare gauge field polarization alters 
only the higher order corrections to the leading S ~ T 2 
scaling of QED3 which itself is determined exclusively by 
the fermion polarization. Equivalently, the finite tem- 
perature modifications to the original QED3 the- 
ory of nodal quasiparticles interacting through vortex- 
antivortex fluctuations are in effect corrections to the UV 
cutoff of the pure 2D quantum electrodynamics. How- 
ever, the leading order scaling of the thermodynamics of 
pure QED3 is dominated entirely by the IR fixed point, 
and is consequently independent of the UV cutoff. 

Spin susceptibility: 
The "topological" fermion spinors tp H allow us to ex- 
press the physical spin density V'jV'T — rt^l as a Dirac 
fermion density xp^fQip. Consequently, the theory con- 
nects correlations among topologically frustrated BdG 
spinons to charge fluctuations in QED 3 . Such correla- 
tions suppress spin response in phase-disordered under- 
doped cuprates and at T = one finds \ ~ Q 2 - Of 
course, within pure QED3 the charge fluctuations would 
still vanish even at finite T since that system is always 
incompressible. Within the QED3 theory of cuprates, 
however, we must now use the "non-relativistic" modi- 
fied finite T bare Lagrangian QJ. As shown bellow, this 
translates into a finite spin susceptibility whose precise 
value reflects the static limit of of the scaling function 
f T . This is in contrast to the thermodynamics where the 
leading low temperature behavior was determined by the 
pure "relativistically invariant" QED3 and is insensitive 
to the specific structure of f T ±. 



To compute the spin-spin correlation function 
(S z (—k)S z (k)) we introduce an auxiliary source J^{k) 
and couple it to fermion three-current. Thus 

C[ip, ip, a M , Jn] = -0 (ry M <9 M + 7p( m + j a0) ^ + £ oW 

(9) 

and since it is the z-component of the spin that couples 
to the gauge field we have 



(S z (-k)S z (k)) 



1 



Z[J^} 6J (-k) SJ a (k) 



Ju=0 



(10) 

Z being the quantum partition function. Now we set 
a' = a M + J M and integrate out both the fermions and 
the gauge field a'. The correlations between a'^ fields are 

) described by the polarization matrix which to the order 
(8) 1/N can be written in the form n M „ = (TI^ + Il^A^ 
(Tin + Ho)B,,,,. The resulting soin correlation functi 



(Ifg + Hg jB^. The resulting spin correlation function 
is then readily found to be 

(S z (-k)S z (k)} = 



ri5 + k' 



(11) 



where 11^ denotes the fermion current polarization func- 
tion. Due to the scale invariance of the (massless) topo- 
logical fermion action, the time component of the re- 
tarded polarization function has the scaling form 



nF ret 
A 



(|k|,o; ) T)=iVTPj(^,- 



(12) 



where P^(x,y) is a universal function of its arguments 
and TV is the number of the four component Dirac fermion 
species (N = 2 for a single Cu02 layer). In the static 
limit, u> — > 0, 



lim lRe P%(x,y) 

3/->0 

and 



lim me P% (x, y) 



21n2 



+ — + 0(x 3 ); x«l (13) 



24tt 



6C(3) 



+ 0(x- 3 ); x>l (14) 



while 



lim 3m P A (x, y) 

y->0 



2 In 2 y , n ty 



(15) 



Furthermore, in order to complete the computation of 
spin susceptibility, we observe that in the limit u> = 
and k — > we have 



fr 



T 



, 00, KT, K T T 



1 



T 2 

'k^ 



(16) 



where c is a pure number [l4j . The expression (|16(l for 
the static limit of f T is rather general and can be un- 
derstood on physical grounds: in the pseudogap state 
at T — > we can think of the "normal" system as a 
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dual superfluid of relativistic "vortex bosons" [5j . These 
bosons are coupled to a massless dual gauge field which 
describes the original long range vortex-vortex interac- 
tions. Once T is finite, this dual superfluid immediately 
losses its off-diagonal order since the screening by the 
dual gauge field allows for creation of a finite density of 
thermally activated dual vortices, i.e. "vortices" in the 
dual superfluid of original vortices. Thus, there is no fi- 
nite T transition in the dual superfluid that would be 
analogous to the Kosterlitz-Thoulcss transition in origi- 
nal superconductor. The above form for f T (| 1 fc>|) is just 
the mathematical expression of this fact. 

Combining all of these results together we finally ob- 
tain the spin susceptibility: 



X (co = 0, |k| 0) 



2Aln2 



K T T 2 



27V In 2 



K T T 



(17) 



Consequently, Wilson ratio R — yT jc v vanishes as T — > 
14J implying the non-Fermi liquid nature of the pseu- 
dogap state in cuprates within the theory of Ref. |5| ■ 

Corrections to scaling: 
There are several sources of corrections to the QED3 scal- 
ing, and it is not practical to address them all within 
the confines of this Letter. Instead, we choose to concen- 
trate here on those particular corrections to scaling which 
arise from the Dirac cone anisotropy. Such corrections 
are most likely to be important in real materials where 
the values of velocities vf and v/\ as functions of under- 
doping are now becoming known with increasing preci- 
sion It was shown in Ref. 0] that the Dirac cone 
anisotropy old = vf/v& scales to unity at the QED3 IR 
fixed point. Alternatively, when defined via an = 1 + 5, 
8 has a small, but finite scaling dimension 775 > 0. To the 
leading order in 1/N, r)s — 32/(5ir 2 N) [J^. For simplic- 
ity, we assume here that the bare action for the gauge 



field is e~ 2 F 2 v 



Then the T— dependent part of the free energy scales as 



(K = e 2 ), where F) 



1"' 



T 3 

T= — -$ N 



T Vp 
2 ' 

jp \C "A/ Lip 

and in the limit of T — > we have 



T 3 fT 
= -^n[-,1 + S(T) 



T 3 

— $jv (0,1) 



— 2 — ^o^jv (0,1) 

Dtp 



(18) 
(19) 



where $at and <fr' N are pure numbers and <5o = an — 1. 
Thus, although the leading order scaling in T is ana- 
lytic, just like in the pure finite-T QED3, the fact that 
the next-to-leading order correction is not reflects the 
strongly interacting nature of the QED 3 IR fixed point. 

In summary, we have shown that the leading order low 
temperature scaling of the thermodynamic functions in 
the pseudogap state is that of a pure finite-T QED 3 , since 
the finite temperature modifications QJ to the theory of 
Ref. |5( affect only corrections to scaling. As a result, the 



specific heat goes as c v ~ T 2 . Furthermore, we argued 
that at low temperatures the spin susceptibility is sup- 
pressed, x i5 T 2 , due to correlations among BdG spinons 
mediated by an emergent gauge field a^. This implies a 
vanishing Wilson ratio and breakdown of the Fermi liq- 
uid behavior within the pseudogap state. Finally, the 
corrections to scaling arising from Dirac cone anisotropy 
were argued to be non-analytic, representing yet another 
manifestation of the strongly interacting character of the 
QED3 brand of quantum criticality. 
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